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Abstract 



Entropy arises in strong interactions by a dynamical separation of "par- 
tons" from unobservable "environment" modes due to confinement. For 
interacting scalar fields we calculate the statistical entropy of the ob- 
servable subsystem. Diagonalizing its density matrix yields field pointer 
states and their probabilities in terms of Wightman functions. It also 
indicates how to calculate a finite geometric entropy proportional to a 
surface area. 
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The long-standing "entropy puzzle" of high-multiplicity events in strong interac- 
tions at high energy has been analysed from a new point of view [I]. The problem 
dates back to Fermi and Landau and is related to understanding the rapid thermal- 
ization of high energy density (^> 1 GeV/fm 3 ) matter |J. Why do thermal models 
work so well? Why do they work at all? 

Or, why does high-energy scattering of pure initial states lend itself to a statis- 
tical description characterized by large apparent entropy from a mixed-state density 
matrix describing intermediate stages in a space-time picture of parton evolution? 
Effectively, unitary time evolution of the observable part of the system breaks down 
in the transition from a quantum mechanically pure initial state to a highly impure 
(more or less thermal) high-multiplicity final state. In Ref. [0] this was discussed 
in detail. Based on analogies with studies of the quantum measurement process 
( "collapse of the wave function" ) |3j and motivated by related problems in quantum 
cosmology and by non- unitary non-equilibrium evolution resulting in string theory [ffl, 
we argued that environment-induced quantum decoherence solves the entropy puzzle 
of strong interactions. 

A complex pure-state quantum system can show quasi-classical behaviour, i.e. an 
impure density (sub)matrix together with decoherence of associated pointer states 
in the observable subsystem |IJ. In particular, there is a Momentum Space Mode 
Separation due to confinement, which is defined in the frame of initial conditions 
for the time evolution and for the physical (gauge) field degrees of freedom. Thus, 
almost constant QCD field configurations form an unobservable environment, since 
they neither hadronize nor initiate hard scattering among themselves. It interacts 
with the observable subsystem composed of partons ||. 

The induced quantum decoherence and entropy production were studied in a 
non-relativistic single-particle model resembling an electron coupled to the quantized 
electromagnetic field, however, with an enhanced oscillator spectral density in the 
infrared. The Feynman- Vernon influence functional technique for quantum Brownian 
motion provided the remarkable result that in the short-time strong- coupling limit 
the model parton behaves like a classical particle ]![]: Gaussian parton wave packets 
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experience friction and localization, i.e. no quantum mechanical spreading, and their 
coherent superpositions decohere. 

Summarizing, partons feel an unobservable (gluonic) environment, which man- 
ifests its strong non-perturbative interactions on a short time scale (-C 1 fm/c) 
through decoherence of partonic pointer states, their quasi- classical behaviour, and 
entropy production. If confirmed in QCD, this will have important consequences for 
parton-model applications to complex hadronic or nuclear reactions 0. The emer- 
gence of structure functions from initial-state wave functions will be further studied 
in our approach. 

We defined a model of two coupled scalar fields representing partons and their non- 
perturbative environment. In the functional Schrodinger picture employing Dirac's 
time-dependent variational principle we derived its Cornwall- Jackiw-Tomboulis (CJT) 
effective action and the equations of motion for renormalizable interactions [1], 
Thus, analysis of the entropy puzzle in strong interactions leads to study an observ- 
able field (open subsystem) interacting with a dynamically hidden one (unobservable 
environment), i.e. quantum field Brownian motion. 

In the following we derive the entropy in any system of two interacting real scalar 
fields. Their most general normalized Gaussian wave functional in the Schrodinger 
picture can be written as 

■ expf-i^-^i^lfGia^-ffii^t)]^-^)]} , (1) 

with (j = 1,2) 

VotWA = ( 2 ) 
Nj(t) exp {-[&■ - fa{t)\ [\Gj\t) - iEj(t)] [fa - fa(t)} + m 3 (t)[fa - fa(t)}} . 

We suppress all spatial integrations. The normalization factors are 

Nj{t) = det{2G,(t)}- 1 / 4 , N 12 (t) = det{l-G 1 (t)G' 12 (t)G' 2 (t)G' 12 (t)} 1 / 4 , (3) 
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discarding an irrelevant constant factor in Nj. Thus, the time-dependent Hartree-Fock 
approximation (TDHF) for the quantum field Schrodinger equation |TJ, [| is embodied 
in the variational parameter one-point functions <f)j(x,t), 7tj(x,t) (mean fields) and 
symmetric two-point functions Gj(x,y,t), Hj(x,y,t), Gi 2 (x,y,t), T,i 2 (x,y,t) (related 
to Wightman functions). Their meaning was discussed in M. 

All physical quantities of the complex system can be calculated with \J>i 2 , express- 
ing inner products by functional integrals. The functional density submatrix p-p for 
the observable "parton" subsystem (0i) is obtained by tracing over the unobservable 
degrees of freedom (0 2 ), 

pp(t) = Tr 2 |*i 2 (t))<* 12 (t)| , (4) 

as calculated explicitly in [|l] (we henceforth omit V). The matrix elements of p 
contain all the information about the subsystem. Our aim is to obtain the von 
Neumann or statistical entropy, S = — Tri pin p. Before, we calculated the simpler 
linear entropy directly, which provides a lower bound for the statistical entropy [0J , 

l-Gi(t)Gi2(t)G 2 (t)Gi 2 (tA ( ^ 
l + G 1 (t)E 12 (t)G 2 m 12 (t)) ' 1 j 

tracing over coordinates. Equation (|5|) is also valid for non-translation invariant 

systems, which is relevant for calculating the geometric entropy related to spatial 

boundaries separating observable and unobservable subsystems. 

Geometric entropy is intimately connected to black-hole entropy ||. Here, one 
identifies <pi as the part of a scalar field <p with support outside a given spatial region 
and 02 = <P ~ 0i) which has its support inside the complement. Our results (|5[) and 
(|T7j) below indicate that geometric entropy comes out finite, once a renormalization 
of the equations for the two-point functions, G's and E's in (g), is performed or a UV 
regularization introduced to provide sufficient integrability constraints. 

We proceed by diagonalizing p. Determining its eigenstates and eigenvalues is 
equivalent to constructing field pointer states [|], 0] within TDHF and their probabil- 
ities. The eigenvalue problem p\p) = p\p) to be solved is of the form 

pF[</>] exp{-0«0 + /30} = (6) 
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S(t) > -- Tr In 



N 2 exp{-0a0 + 60} J V<f>' F[0'] exp{-0'[a* + a]0' + [b* + (3 + 0c*]0'} , 

using the ansatz ((0 + 0i)|p) = F[(f)] exp{— 0a0 + /3(f)} with unknown one- and two- 
point functions (3 and a and a non-exponential functional F . According to results for 
p-p from I], we define iV = N\Ni 2 , b = mi (0 X does not appear in (||)), 

a = iGT 1 A-2[S 1 -|(G' 12 G' 2 S 12 + S 12 G' 2 G' 12 )] = a* , 
c ee 1 -B — |[Gi 2 G 2 S 12 — S 12 G 2 Gi 2 ] = , 

and combinations of two-point functions 

A ee 1 — |GiGi 2 G 2 Gi 2 + ^GiSi 2 G 2 Si 2 , 
-B ee ^GiGi 2 G 2 Gi 2 + |GiSi 2 G 2 Si 2 . 

Choosing (3 = b in (^), completing the square, shifting 0', and requiring resulting 
Gaussians in to cancel yields the eigenvalue problem: 

p F[<j>] = N 2 J V<f>' F[<f>' + Y<j>\ exp{-0'X0'} , (7) 

with X = a* + a = X 1 , Y = |X _1 c, and where a = a 1 , by (|6]), is determined to solve 
the equation a — a = |c*[a* + a] _1 c. Note the similarity to the finite-dimensional 
oscillator problem of Srednicki || . 

Equivalents, replacing F[<f/\ -> F[8/8{(j)c*) + 6/6{af>)] and 0c* 0' -> i[0c*0' + 0'c0] 
in (^), we obtain by integration 

p F[0] = N 2 detiX}- 1 / 2 exp{-i0c*X- 1 c0} + j^} exp{+i0c*X- 1 c0} , 

(8) 

which is more convenient than (|7]). Looking for polynomial functional solutions of 
(||), we find first of all a constant, 

F o [0] = 1 po = N 2 detiX}- 1 / 2 . (9) 

Secondly, instead of a general linear functional, the Fourier transform is sufficient, 

F k [<j>] = J d d x e- lkx <P(x) ee 0, , (10) 
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since the problem is linear in F. Then, from (H) — ([10]) 



Pi^k = jPo[X 1 c<p + <pc*X % = p Y kk ,<i)_ k , , (11) 



summing over indices occurring twice. For a translation-invariant system, ( |TT|) could 
immediately be solved. Generally, however, denoting eigenvalues and eigenvectors 
of (Y kk r) by £fc and <p k , one obtains a set of linear eigenvalues p k = po6c- Due to 
the Gaussian structure in (H), the higher-order eigenfunctionals can be built up as 
linear combinations of products of 4> k s and lower-order ones. For example, F kk /[(f)] = 
4>k4>k' + C kk >, which yields a set of quadratic eigenvalues p kk * = Po6c£fc'@(^' ~ k). Note 
the constraint k! > k; interchange of k and k! does not lead to a new eigenfunctional 
due to the scalar (bosonic) character of the fields. The constant C kk i follows with the 
help of the matrix diagonalizing (Y kk >). We do not construct explicitly the higher- 
order eigenfunctionals. However, the n-th order set of eigenvalues, 

n 

Pki...k n = Mfc! II £fci©(^ - k^) , (12) 

i=2 



is easily found, similarly to p kk > above. To check the result (|12|), we calculate 



Tr = J2 eigenvalues = p + J2 Pki-k n 

n=l ki...k n 

= po[i + £*i + Yl &i&a + Yl + • • • 

fcl fci<fc2 ki<k2<k'j 
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= Poll E 4"* = Pont 1 -^]" 1 = Podetjl-rr 1 = 1 , (13) 

k rifc=0 k 

which resembles the evaluation of a bosonic partition function. In the last step we 
used po = det{X- 1 Re[2a - c}} 1 ' 2 = det{[l - \X- X (f][\ - ^X^c}} 1 / 2 , which follows 
from the equation determining a or X. 
Similarly, we obtain the linear entropy, 

S Un = Tnfp-p 2 } = 1- T ri p 2 = 1- det{[^} . (14) 

In order to express Y in terms of A and B, we observe that in a direct calculation Q 
of Trip 2 (and in the n-fold functional integral for Trip n ) imaginary parts of a and c 
cancel. Therefore, we replace a and c by their real parts, a = jG^A and c = ^G^B, 



simplifying the equation for a, X, or Y, a — a = \c[a + a] 1 c. The solution (for 
integrable eigenfunctionals) is 

Y = ( C /2)- 1 /2[ a + (a 2_ l) l/2j- 1 (c/2) l/2 ; (15) 

with a = (c/2y 1/2 a(c/2y 1/2 . Finally, inserting (ITJ) into (p 



which confirms our earlier result, employed in (HD Q]. Next, we calculate the statistical 
entropy using the "replica trick" : 

= -Tr|ln(l-y) + ^— ylnyj . (17) 

Together with flTSP, eq. (0) presents our main result. It generalizes eq. (^) of 
Srednicki ||. Basically, the TDHF approximation for interacting quantum fields 
preserves a Gaussian structure of the wave functionals, see (1) — (3), which is exact in 
the non-interacting case and can be reduced to a coupled harmonic oscillator problem. 

To evaluate the entropy ( |17|) is still a formidable task for any realistic situation. 
Before trying, it seems worth while to draw some general conclusions: 

I. Neither mean fields 0^2, nor their conjugate momenta 7^2, nor imaginary parts 
Si_2 of the "parton" and environment two-point functions contribute to S. 

II. Vanishing correlations between "partons" and environment, i.e. G12 = £12 = 
(independent subsystems), imply A = 1, B = 0, i.e. Y = 0, and S = 0. 

III. Vanishing widths of "parton" or environment wave functionals, i.e. — > 
(one or the other subsystem classical/reversible |I]), imply Y = and S = 0. 

This presumably holds for any field theory of "partons" coupled to environment modes 
independently of the interactions in TDHF approximation. The time-evolution, how- 
ever, follows specific equations of motion for the one- and two-point functions JT|. 

Our considerations confirm that quantum decoherence and entropy production in 
a subsystem is induced by an active environment ]IJ [|, |, || . The above diagonaliza- 
tion of the "parton" density functional yields time-dependent field pointer states, the 



simplest one of largest probability po being 



# [<M = exp{-[0-0i(t)]a(*)[0-^i(t)]+i7fi(*)[0-^i(*)]} , (18) 
cf. (H) — with a = [(c/2)(a 2 — l)(c/2)] 1 / 2 . Higher-order eigenfunctionals are less 



probable, see (12), and have higher kinetic energy, since their wave functionals have 
additional nodes, e.g. fllOD , analogous to excited oscillator states. 

As a first application of ( fL7j) we consider the large-entropy limit, i.e. Y « 1 or 
A « 5. Then, using ( |i4|) and (|T6|), we find: 

S(t) « -Tr ln(l-y)/(l + y) = -±Tr \n(A - B) / (A + B) 

— |E^=i ^ Tr {[GiGi2G2Gi2] n — [— GiSi2G2Si2] n } , (19) 

i.e. (0). If we assume a spatial surface of area .4 dividing the system into two, 
which is flat on the scale of the short-ranged correlations in (|19"D , then Tr[. . can be 
interpreted as a sum of closed loops of strings of G"s or S's intersecting the surface 2n 
times: once for each factor G 12 or S 12 correlating in- and outside fields. The dominant 
contribution to the trace comes from small loops (let <C G± 2 )- Regularizing their 
contribution by a short- distance cut-off C, their size transverse to the surface will be 
0(C 2 ) for d = 3. Transverse to the surface the system is locally translation- invariant. 
Therefore, the geometric entropy is approximately 

where Tr^ is evaluated locally on the scale of C. A dimensional analysis led Srednicki 
to propose S oc A before ||. Equation fl2~0|) can be applied to the moving mirror 
model; following Kabat et al. 0, it approximates the thermal entropy outside a black 
hole of radius much larger than C 

Secondly, coming back to partons interacting with their (gluonic) environment, 
the rate of entropy production, which follows from (JPTD, is most interesting. We define 
a dynamical decoherence time r, 

W dt W Tr yiny JdkY k \nY k V ; 



with dk = d d k/(2ir) d . For simplicity we assumed small Y or S and a translation- 
invariant system; the Fourier transform is Y\. = B^[A^ + (A 2 — B^) 1 ' 2 ]" 1 , since A, B 
are convolutions of two-point functions now. Generally, two limits are particularly 
important: r(t — > 0) gives the time scale for the decay of a Gaussian partonic field 
state, cf. (0) — (|3|), into an incoherent superposition of pointer states, e.g. fllgp, with 
impure density matrix and non-zero entropy; r(t ^> 0) reflects the approach to a 
stationary state (thermalization) , if it exists. Using the equations of motion from 
the decoherence time will be calculated for phenomenologically interesting situations 
elsewhere. 

I thank N. E. Mavromatos, L. Pesce, and J. Rafelski for stimulating discussions. 
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